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Abstract 

The hyperfine splitting in hydrogenlike 209 Bi, 203 T1, and 205 T1 is calculated with the nuclear 
magnetization determined from experimental data on the hyperfine splitting in the corresponding 
muonic atoms. The single-particle and configuration-mixing nuclear models are considered. The 
QED corrections are taken into account for both electronic and muonic atoms. The obtained re- 
sults are compared with other calculations and with experiment. 
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1 Introduction 



High-precision measurements of the hyperfine splitting (HFS) in heavy hydrogenlike ions fflEl 
Ell31|5l have triggered a great interest to calculations of this effect. The main goal of these exper- 
iments was to probe the magnetic sector of quantum electrodynamics (QED) in the presence of a 
strong Coulomb field. The uncertainty of the theoretical results is mainly determined by the un- 
certainty of the nuclear magnetization distribution correction, the so-called Bohr-Weisskopf (BW) 
effect. In calculations, based on the single-particle nuclear model lIUIZllIDE, w ^i c h provide a rea- 
sonable agreement with the experiments, this uncertainty may amount up to about 100% of the 
BW effect and is generally few times larger than the total QED contribution. More elaborated cal- 
culations, based on many-particle nuclear models ifTOl ITT1 . do not provide a desirable agreement 
with the experiments. 

In the present paper, we determine the BW correction to the hyperfine splitting of hydrogen- 
like 209 Bi, 203 T1, and 205 T1 using experimental data on the hyperfine splitting in the corresponding 
muonic atoms. We consider the single-particle and configuration-mixing nuclear models. The pa- 
rameters of the nuclear magnetization distribution are chosen to reproduce the experimental values 
of the nuclear magnetic moment as well as the BW contribution in muonic atoms extracted from 
the corresponding experiments. To increase the precision of determining the BW contribution, the 
QED corrections for the muonic atoms have been evaluated. The obtained results are compared 
with other calculations and with experiment. 

The relativistic units (h = c = 1) and the Heaviside charge unit (a = e 2 /(47r), e < 0) are 
used in the paper. 
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2 Hyperhne splitting in muomc atoms 



The ground-state hyperfine splitting in muonic atoms can be written in the form: 

AE = AEns + AEbw + AEqed, (2.1) 

where AE^s is the hyperfine splitting value incorporating the relativistic and nuclear charge dis- 
tribution ("nuclear size") effects, AEbw is the BW contribution, AEqed is the QED correction. 
The AE^s value can be calculated by the formula: 

oo 

AE m = - a iJL±!E+ll [ dr g(r)f(r), (2.2) 
3 fj, N m p 21 J 
o 

where a is the fine structure constant, fi is the nuclear magnetic moment, /x/y is the nuclear 
magneton, m p is the proton mass, and I is the nuclear spin. g(r) and f(r) are the radial parts of 
the Dirac wave function: 

which are determined by solving the Dirac equation with the Fermi distribution of the nuclear 
charge (47r f drr 2 p{r) = 1): 

P(r) = — P0 (r . cy (2.4) 
1 + exp(— ) 

Here c is the half-density radius and a is related to the skin thickness t by t = (4 log 3)a, defined 
as the distance over which the charge density falls from 90% to 10% of its maximum value. 

The individual contributions to AE for muonic atoms of 203 T1, 205 T1, and 209 Bi are presented 
in Tabled The AE^s values are given in the second column. In the third column we present the 
BW correction evaluated within the single-particle nuclear model according to the prescriptions 
given in |6] El HI (see also the next section of the present paper). The wave function of the odd 
nucleon was determined by solving the Schrodinger equation with the Woods-Saxon potential. In 
muonic atoms, the QED correction is mainly determined by the vacuum polarization (VP) con- 
tribution, which consists of the electronic electric-loop and magnetic-loop parts. In the Uehling 
approximation, the electric-loop part is determined by the potential: 

oo oo 

CC(r) = -aZ 2 ^ J dr'^p{r') f dt (l + -L) ^El (2.5) 

o 1 
exp(— 2m\r — r'\t) — exp(— 2m{r + r')t) 
4mrt 

where m is the electron mass. The corresponding correction (AEyp) is derived as the difference 
of equations (12.21) with the wave functions obtained by solving the Dirac equation with and without 
the Uehling potential J2.5I) . For the correction to the hyperfine splitting due to the magnetic loop 
one obtains: 

AE^ = (A\U^(r)\A), (2.6) 
where \A) is the state vector of the whole atomic system, 

oo 

U^(r) = ^hfs(r)^ Jdt(l + Ej ^ t2 ~ 1 (1 + 2mrt) exp(-2mrt), (2.7) 
l 

i?hfs(r) is the hyperfine interaction operator: 

_ je Ka-fr xr]) 
47r r 6 
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Table 1: Individual contributions to the hyperfine splitting in muonic atoms, in keV. 



Atom 












203 


4.70 


-2.06(66) 


0.06 


2.70(66) 


2.66(30) L12J 


205 


4.73 


-2.06(66) 


0.06 


2.73(66) 


2.32(6) Q 


209 Bi 


6.69 


-2.49(80) 


0.09 


4.29(80) 


4.44(15) O 



H is the nuclear-magnetic-moment operator, and a is a vector incorporating the Dirac matrices. 
The total QED correction (AEqed ~ ^-Eyp + ) is § iven in tne fourt h column of Ta- 

ble [2 The total theoretical values obtained in this work (fifth column) are in a good agreement 
with the experimental ones (sixth column) and with the previous theoretical calculations [14|, 
which do not account for the QED corrections. Since the theoretical uncertainty is mainly deter- 
mined by the uncertainty of the BW effect, the experimental values of the hyperfine splitting in 
muonic atoms can be employed to determine the BW contribution and, therefore, the parameters 
of the nuclear magnetization distribution for a given nuclear model. Then, with these parameters, 
the BW correction to the hyperfine splitting in the corresponding electronic ions can be calcu- 
lated. Such calculations are presented in the third and fourth sections for the single-particle and 
configuration-mixing nuclear models, respectively. 



3 Nuclear magnetization in the single-particle model 



In the single particle model, the nuclear magnetization is ascribed to the odd nucleon. The nuclear 
magnetic moment is given by 



li \\[g s + (2I-l)g L ], 



I = L + l/2 



UN 



t— l-g s + (21 + 3)g L ] , I = L- 1/2, 



(3.1) 



2(7+1) 



where I and L are the total and orbital angular momenta of the odd nucleon. For proton gi = 1 
and for neutron g^ = 0. In calculations of the HFS the g$ factor is usually chosen to yield the 
observed value of the nuclear magnetic moment. To calculate the BW effect within the single- 
particle model, one has to adopt the replacement n — > n(r) = F(r)n in the HFS operator. The 
function F[r) is given by 151 : 



F(r) 



fJ>N f 



\ 9 s + (i-\Vjl 



ar r u (r ) 



+ 



21-1 ( 1 

sTTTl) 55 + V 2 ' ' //; 



drW{r>) (^) 3 } 



(3.2) 



for / = L + 1/2 and 

Fir) ■ 



Mat f 
H I 



2(1 + 1 



/ 1(21 + 3) 

9S + „, r , J ql 



2(1 + 1) 



j / 12 2/ l\ 

ar r u [r ) 



+ 



21 + 3 1(21 + 3) 

:9S + K7T , TT 9L 



8(1 + 1) 



2(7 + 1) 



oo 

j dr'r' 2 u 2 (r') (^) 3 } (3.3) 
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Table 2: The Bohr-Weisskopf correction to the HFS in electronic hydrogenlike ions, derived from 
experiments on muonic atoms within the single-particle nuclear model. fe™° d denotes the uncertainty 
of e s p caused by the different parameterizations used in the calculation (see equations (13.51) . (I3.6!) ). 



Ion 


r 


Z^ xp (keV) 


^s.p. 


r_mod. 
ofc s.p. 


203^80+ 


5+ 


2.66(30) [12 1 


0.0155(35) 


9% 


205rpi80+ 


1+ 


2.32(6) d 


0.0193(24) 


11% 


209g|82+ 


9 
2 


4.44(15) O 


0.0123(14) 


8% 



for / = L — 1/2. Here u(r) is the radial part of the wave function of the odd nucleon. The relative 
value of the BW correction, defined by e = — Z\-Ebw/^-^ns> can be written in the form: 

oo 

J drF{r)g{r)f(r) 

£s. P . = 1 - ^ • (3-4) 

/ drg(r)f(r) 
o 

To reproduce the experimental values for the HFS in muonic atoms within the single-particle 
nuclear model, we have considered the following parameterization for u(r): 

( u (r) = u r n , n = 0,1,2, r < R M „ 
\u(r) = 0, r > R M 

and 

fu(r) = u Q (R M - r) n , n = 1,2, r < R M 
\u(r) = 0, r>R M . 

The constant no is determined by: f °° drr 2 u 2 (r) = 1. The magnetic radius i?M is derived from 
the equation: 

Z\£ c % - AE QED = (1 - £ s . p .(#m))^£ns, (3.7) 

where AEe xp and Z\£qed are the experimental value of the HFS for muonic atom and the QED 
correction, accordingly. Then we can calculate the BW correction for the corresponding electronic 
ions with Rm derived from equation i3.lt . We have found that the results for the BW correction 
(e s . p .) in electronic H-like ions are stable enough within the parameterizations (13.5b . (I3.6I) . For this 
reason, the uncertainties of e s . p , presented in Table|2j are mainly determined by the uncertainty of 
the experimental HFS values in muonic atoms. 



4 Nuclear magnetization in the configuration-mixing model 

In the configuration-mixing model the nuclear magnetism is determined by the last odd nucleon 
and the particle-hole excited states. So, the HFS (without the QED correction) can be represented 
as: 

AE c . m . = AE s . p .+5E{Afi), (4.1) 

where A\i is the correction to the nuclear magnetic moment, due to mixing particle-hole states. 
The nuclear magnetic moment can be written in the form: 

«cx P = «s. P . + Aii. (4.2) 
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The formulas for AE S . P . are well known (see, e.g., [ 14]): 

4 1 21 + 1 



AE< 



s.p. 



-a- 



3 m p 21 



\gs+(l-\)c lL 



s.p. 



-a 



+ 



4 1 21+1 
3 m p 21 



21-1 / 1 

L WTT) 9S + V-2 },JL 



K a + 

K b ], I = L + l/2, 



(4.3) 



/ 1(21 + 3) 

rffs + -^7T—rr9L 



+ 



2(1 + 1)*° 2(1 + 1) 
21 + 3 1(21 + 3) 



:9s + 



2(1 + 1: 



-9L 



where 



L8(/ + l 

oo oo 

K a = I dRR 2 ul L (R) J drg(r)f(r), 



K b ], I = L-l/2, 



o 

oo 



R 



R 



K b = J dR R~ l u 2 n L (R) J drg(r)f(r). 



(4.4) 



(4.5) 



(4.6) 



The correction for the AL = mixing terms is given by fiM : 
5E(Afi) = 



at 1 2I „t 1 / dr f(r)g(r)^2( L ,A(i 

L' 



3 m p 21 



o 



l 



{l-K s } Ll +^-^{K s -K L } L , 



9S ~ 9L 



l' 



(4.7) 



(4.8) 



The quantities C, L > determine the re-distribution Afi between the configurations. We use the fol- 
lowing designations: 



oo 

{A} L , = J dRR 2 u n)L/;J (R)u n L , )f (R)A(R), 



(4.9) 



K S (R) 



J drg(r)f(r) 
b 

oo 

/ drg(r)f(r) 



(4.10) 



R 



fdr(l-^)g(r)f( 



K L (R) 



J drg(r)f(r) 




(4.11) 



Here u n ^ and it n ,L,j are the radial parts of the wave functions of the odd nucleon. They are 
obtained by solving the Schrodinger equation with the Woods-Saxon potential. We consider Q L i = 
jj, where M is the number of the configurations. The quantities Z\/i and g$ are derived from the 
equations: 



^ex P _ )\9s + {I-\)9l + ^T.l'Cl'Wl^ I = L + l/2 

WW)9s + ff^9L + ^E L 'C L '{l} L ', I = L-l/2, 



UN 



(4.12) 
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Table 3: The Bohr-Weisskopf correction to the HFS in electronic hydrogenlike ions, derived from 
experiments on muonic atoms within the configuration-mixing model. 



Ion 


r 




^c.m. 


203ji80+ 


\+ 


2.66(30) L12J 


0.0179(36) 


205-Q80+ 


1+ 


2.32(6) (El 


0.0214(6) 


209g|82+ 


9 
2 


4.44(15) [13| 


0.0119(11) 



AEgp - AE QED = AE S . P . (g s ) + 5E {Api) . (4.13) 
The Bohr-Weisskopf correction is determined by: 

_ AE s . P X9s) + SE(Afi) 

Ec - m -~ AE m ■ (4 - i4) 



With gs and A^i derived from equations (I4.12t . (I4.13t . we can calculate the BW correction for 
the electronic ions. The results of these calculations are presented in Table|3] The values of e c m . 
given in Tableware in a good agreement with e s . p . presented in Table|2] 



5 Results and discussion 



In Table |4] we compare the BW correction to the HFS in electronic H-like ions, derived from 
the experiments on muonic atoms. e s p . and e c . m . are obtained employing the single-particle and 
configuration-mixing models, respectively, e is a value, which corresponds to the parameterization 
d3.5t with n = and whose uncertainty covers all the e s p . and e c . m . values. For comparison, the 
values of the BW correction obtained previously by direct calculations within the single-particle 
and many-particle models are presented as well. We conclude that our results for e are stable 
enough with respect to a change of the nuclear model. They also have a better accuracy than the 
previous single-particle results 1^1 d [8] |9j . 

In Table[5]we present our final theoretical results for the HFS in electronic hydrogen-like ions. 
These results, which include the BW correction derived in the present work and the QED correc- 
tion taken from [TJIT^inilEL wt compared with previous calculations and with experiment. As 
one can see from the table, our results are closer to the experimental ones, compared to the results 
based on the direct calculations within the single-particle nuclear model [8|. However, due to a 
higher accuracy of the present results, a small discrepancy between the theory and experiment 
occurs for 209 Bi 82 +. The reason for this discrepancy is unclear to us. 
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Table 4: The Bohr-Weisskopf correction to the HFS in electronic H-like ions, derived from experi- 
ments on muonic atoms. 









209gj82+ 


205rpj80+ 








^s.p. 


0.0123(14) 


0.0193(24) 


0.0155(35) 






^x.m. 


0.0119(11) 


0.0214(6) 


0.0179(36) 






e 


0.0123(15) 


0.0193(27) 


0.0155(40) 


Shabaev et al. 




e 


0.0118 


0.0179 


0.0179 


Labzowsky et al. 





e 


0.0131 






Tomaselli et al. 


Co| 


e 


0.0210 






Sen'kov and Dmitriev 


LI U 


£ 


0.0095 (+ 7 8 ) 







Table 5: The total theoretical results for the hyperfine splitting in electronic H-like ions, in eV. 



Theory Theory Experiment 

[this work] [8| 

TP+ 3.220(20) 3.229(17) 3.21351(25) 

Tl 80+ 3.238(9) 3.261(18) 3.24410(29)0 

Bi 82+ 5.098(7) 5.101(27) 5.0840(8) [1] 
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